Complex Solitary Waves and Soliton Trains in KdV and mKdV Equations 
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We demonstrate the existence of complex solitary wave and periodic solutions of the Kortweg de- 
vries (KdV) and modified Kortweg de-Vries (mKdV) equations. The solutions of the KdV (mKdV) 
equation appear in complex-conjugate pairs and are even (odd) under the simultaneous actions of 
parity (V) and time-reversal (T) operations. The corresponding localized solitons are hydrodynamic 
analogs of Bloch soliton in magnetic system, with asymptotically vanishing intensity. The PT-odd 
complex soliton solution is shown to be iso-spectrally connected to the fundamental sech^ solution 
through supersymmetry. 


I. INTRODUCTION 

The celebrated Kortweg-de Vries (KdV) equation is 
a well-studied non-linear dynamical system, first evoked 
for the description of solitary waves [T] in shallow water 
[ 2 ]. It has, since then, found much applications [3]. An 
integrable model [4] , it arises from the compatibility con¬ 
dition of two linear equations, the well-known Lax-pair 
[5]. The fact that one of them is the linear Schrodinger 
operator, connects this non-linear equation to quantum 
mechanical eigenvalue problem of the reflectionless po¬ 
tential Ii[7]. It has both localized and periodic cnoidal 
wave solutions, which appear in the Lax equation as po¬ 
tentials, giving rise to bound states [ 8 ] and band structure 
[9], respectively. A number of methods, viz., inverse scat¬ 
tering m, Hirota bilinear m etc., have been developed 
to find the general multi-soliton solutions of the KdV hi¬ 
erarchy. The KdV solutions are connected by the Miura 
transformation m to the modified KdV (mKdV) equa¬ 
tion m, which also has found diverse physical applica¬ 
tions. It appears in the description of van Alfven waves in 
collisionless plasma El, phonons in anharmonic lattice 
m, interfacial waves in two-layer liquid with gradually 
varying depth m, transmission lines in Schottky barrier 
m and ion acoustic solitons [T8]-[20], to mention a few. 

In addition to the fundamental localized soliton so¬ 
lution: —2a^sech‘^[a{x — da^t)] and its periodic coun¬ 
terpart, —21710^cn^[a{x — 4(2m — l)a‘^t)], recently real 
superposed solutions have been found to satisfy the KdV 
equation EU Eg. Here, we show that KdV equation 
possesses T^T-symmetric, localized and cnoidal wave so¬ 
lutions, which appear in complex conjugate pairs. Inter¬ 
estingly, the sum of the paired solution is also a solution 
of the KdV equation, whereas the difference is not. The 
existence of these solutions can be traced to the fact that, 
there are two distinct mKdV equations, 
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whose solutions are related through v ^ iv. The Miura 
transformation for the solution of KdV, u = ± 

then implies that, u = —v^ ± ivx, is also a solution. We 
find complex VT-odd solutions for the mKdV equation, 
where the sum and the differences of the pair are also 
solutions. This, in turn, generates more general complex 
superposed solutions for the KdV equation. The complex 
soliton solutions are analogs of Bloch solitons in mag¬ 
netic systems [23]. The T^T-symmetric soliton solution 
is shown to be isospectrally related to the fundamental 
solution sech^x, in the Lax equation. 


II. COMPLEX SUPERPOSED SOLUTIONS 

We start here with the KdV equation, 
duUx T Uxxx — d? 

where, u = u{x,t), = ||, Wa, = and 


A. cn^ ± isndn type solution 

It is straightforward to check that, the following pair 
of complex periodic solutions satisfy the KdV equation: 

u{x, t) = Acn‘^{(, m) + iBsn{(, m)dn{(, m), 

where ( = a{x — ca^t), provided A = —ma^, B = 
2 z^/ma‘^ and c = (2m — 1), m being the modulus param¬ 
eter. Evidently, velocity c exhibits two disjoint domains: 
for ^ < m < 1, the solution is right moving, while for 
the remaining half, 0 < m < it is left moving. It is 
interesting to note that, the sum of the complex pair is 
the well-known cnoidal wave solution, with c = 4(2m — 1) 
and A = —2ma^, whereas the difference of the pair is no 
longer a solution. The real (ur) and imaginary parts (uj) 
of the solutions are plotted in Fig. 1, which are even and 
odd functions of the argument, respectively. The funda¬ 
mental soliton has higher intensity in comparison to the 
superposed solution. 
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FIG. 1: Plots of real (solid) and imaginary (dashed) parts 
of u{Q = cin?{C,) + isn{Qdn{Q ^ for m = 1 (purple) and 
m = 0.25 (green). The inset shows intensities of su¬ 
perposed (red) and fundamental (blue) soliton solutions, 
with the latter being larger. For all the cases, a = 1. 

B. c'n? di isncn type solutions 

The following factorizable, superposed solution also 
satisfies KdV equation: 

u{x, t) = m) + iBsn{(^^ m)cn{(^, m) + 

provided A = —ma^^ B = ±^4 and c = (5m —4) —6/3. As 
in the previous case, here also, oppositely propagating 
modes occupy two different domains of m. The real 
and imaginary parts of the solutions are plotted in Fig. 
2. Here also, addition of the paired localized solutions 
(m = 1) yields the fundamental soliton with higher 
intensity. 
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FIG. 2: Plots of real (solid) and imaginary (dashed) 
parts of u{() = cn‘^{C) + isn{()cn{()^ for m = 1 (purple) 
and m = 0.25 (green). The real parts are symmetric, 
whereas the imaginary ones are antisymmetric about 


For Of = 1, the inset shows intensities of superposed (red) 
and fundamental (blue) soliton solutions, with the latter 
being larger. 

C. Remarks and observations 

For m = 1, superposed pairs reduce to T^T-symmetric 
complex form: u{x;t) = sech^C ± isechCtanhC. Use of 
the Gole-Hopf transformation: = ^, in one of the 

Miura route: u = v‘^ Vx^ yields u = The Galilean 
invariance of the KdV equation [24] allows a constant 
shift in u, leading to 

-^xx + [\^u{xA)]i^ = 0 , 

the one dimensional Schrodinger equation with a Scarf- 
type T^T-symmetric potential [25ll27] . Soliton solutions 
can be generated through iso-spectral deformation of 
the potential [2H1 EH], wherein both the wave function 
and the potential can change, leaving the spectrum 
invariant. Depending on the number of bound states, 
this iso-spectral flow introduces parameters, suitably 
interpretable as time variables m- The complex Scarf 
potential is known to be iso-spectral to the real sech‘^x 
potential [31], explaining the common real eigenvalues 
for both these potentials. 

KdV also posseses real superposed singular [32] so¬ 
lution: u{xA) = A cosech^C + HcosechCcoth^, with 
A = B = and c = 1. 

III. COMPLEX SUPERPOSITION SOLUTIONS 

For the mKdV equation, 

Vt - Qv'^Vx + Vxxx = 0, 

there exist solutions in the form of complex superposi¬ 
tion. 


A. sn ± icn type solutions 

The following parity odd superposition solution solves 
mKdV equation: 

V = Aasn{(, m) + iBacn{(, m), 

provided A = B = a.nd c = f - 1. Unlike 

KdV, the sum of the pair, as well as their difference, 
simultaneously satisfy the mKdV dynamics. In case of 
the sum, v = 2 Aasn{C,^ m) is an exact solution, when 
c = 5(m — 5) and A = while for the difference, 

V = 2iBacn((,m) is also an exact solution, provided 
c = (2m—1) and B = These form of solutions gen¬ 

erate more general superposed solutions for KdV equa¬ 
tion. The real and imaginary parts of the paired solutions 
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are plotted in Fig. 3, which are odd and even functions of 

the argument, respectively. The corresponding intensity 
2 

is constant, / = 



FIG. 3: Plots of real (solid) and imaginary (dashed) parts 
of u{Q = sn{()-\-icn{()^ for m = 1 (purple) and m = 0.25 
(green). The real parts are antisymmetric, whereas the 
imaginary ones are symmetric about The inset shows 
intensities of individual (black), added (red) and sub¬ 
tracted (blue) fundamental soliton solutions, with a = 1. 

B. sn ± idn type solutions 

Another form of complex periodic pair will satisfy the 
mKdV equation: 

V = Aasn{(, m) + iBadn{(^ m), 

provided A = and c = ^ — m. In case of 

m = 1, one obtains v{xA) = ^{tanhC ±isech(), which is 

odd under VT operation. It is interesting to note that, 

the velocity of the pair, sn^icn is —1/2 times the velocity 

of dn solution and the velocity of the pair, sn ± idn is 

— 1/2 times the velocity of cn solution. Fig. 4 plots the 

real and imaginary parts of the paired solutions, which 

again are odd and even functions, respectively. Intensity 

2 

associated with the excitation is constant, / = ^ and 
is independent of the modulus parameter. Real, singular 
superposition solutions of the form v{xA) = A cosechC + 
B cothC, also satisfy mKdV dynamics, with B = A = 
and c = —1/2. 
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FIG. 4: Plots of real (solid) and imaginary (dashed) parts 
of 14(C) = sn{Q^idn{()^ for m = 1 (purple) andm = 0.25 
(green). The real parts are antisymmetric, whereas the 
imaginary ones are symmetric about The inset shows 
intensities of individual (black), added (red) and sub¬ 
tracted (blue) fundamental soliton solutions, with a = 1 . 

IV. CONCLUSION 

In conclusion, trivially T^T-symmetric KdV equation is 
shown to posses T^T-symmetric complex solutions, with 
asymptotically vanishing intensity for the solitons. In 
the corresponding Lax equation, the T^T-symmetric po¬ 
tential supports real eigenvalues and is iso-spectrally con¬ 
nected to the real potential sech^x, the fundamental soli- 
ton solution. Iso-spectral deformation [29] has been use¬ 
ful in generating the multi-soliton solution for the KdV 
equation. The realization of the same in the case of 
complex T^T-symmetric potentials needs careful investi¬ 
gation. The fact that in the T^T-symmetric phase , an 
inner product [33] is defined, may facilitate in obtain¬ 
ing these multi-soliton solutions through iso-spectral de¬ 
formation. For the mKdV equation, VT-odd solutions 
have been found to be exact solutions, which generate 
general T^T-symmetric potentials for the KdV equation, 
through Miura transformation. The KdV and other in- 
tegrable systems, like Boussinesq hierarchy, manifest in 
two-dimensional induced gravity j34ll36] and conformal 
field theory m- The implication of the complex VT- 
symmetric solutions in the context of two-dimensional 
gravity is worth exploring. 
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